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Abstract Under fixed cutting conditions, the surface finish 
roughness is correlated to the grain angle. However, the 
means of determining the grain angle automatically and 
accurately is still a challenge for on-line control of the 
router. It is therefore necessary to develop a new technol¬ 
ogy to determine the grain angle accurately and automati¬ 
cally. In this research, a laser light scattering pattern was 
used to accurately determine the grain angle. The light scat¬ 
tering pattern image was a quasi-ellipse caused by the grain 
direction and tracheid effect. A new modified Hough trans¬ 
form ellipse analysis technology was adopted to determine 
the ellipse parameters that could be used to determine the 
grain angle. The results indicated that the measured grain 
angle using the method proposed here was accurate and 
effective. The measured grain angle coincided with the real 
grain angle. There was an insignificant difference between 
the measured grain angle of Japanese beech (Fagus crenata 
blume) and that of sugi (Cryptomeria japonica D.Don) 
under two machining conditions that gave planed or sawn 
finishes. However, the accuracy of the measured grain angle 
of sugi was better than that of Japanese beech for the 
planed finish, the accuracy of the measured grain angle of 
Japanese beech was better than that of sugi for the sawn 
finish, and the accuracy of the measured grain angle for 
planed samples was better than that for sawn samples of 
both sugi and Japanese beech. 

Key words Grain angle • Laser light scattering pattern • 
Image processing • Ellipse analysis 


Introduction 

Due to increasing labor and raw material costs and high 
customer expectations, it is essential for sawmills to con- 
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tinually improve productivity and quality. The simplest way 
to achieve this goal is to increase the cutting speed. How¬ 
ever, excessive cutting speed will generate a coarse surface. 
Cutting force and acoustic emission (AE) have been used to 
evaluate the grain angle used to control the feed rate to 
achieve the maximum cutting efficiency and finished surface 
quality in an on-line controlled router. 1,2 However, the cut¬ 
ting force is affected by its additional inertia and the differ¬ 
ence between the latewood and earlywood. In addition, the 
linear relationship between AE count rate and surface 
roughness under every discrete grain angle is complex and a 
means of determining the grain angle automatically and 
accurately is still not clear. It is therefore necessary to de¬ 
velop a new technology to determine the grain angle accu¬ 
rately and automatically in the control system of the router. 

Determining the grain angle is a valuable but difficult 
job. Methods based on the dielectric constant and ultrasonic 
wave propagating model for across and along the grain have 
been developed. 3,4 Furthermore, wood cells (either the tra¬ 
cheid cells in softwoods or the fibrous cells in hardwoods) 
were modeled as a bundle of end-to-end tubes lying parallel 
to the stem direction. When the laser beam hits wood sur¬ 
face, it will be transmitted inside the wood materials and 
reflected from the wood surface. In the process of propagat¬ 
ing inside, some of the transmitted light will emerge from 
the wood surface in the area surrounding the incident point 
and form a quasi-ellipse with its major axis along the grain 
direction. The phenomenon has been called the tracheid 
effect because it was first observed in the tracheid structures 
and is more obvious in softwoods than in hardwoods. A 
laser scanning system was developed for wood defects and 
grain direction detection using the tracheid effect. In this 
system, a rotating mirror drum was used to deflect the laser 
beam onto the wood surface and a photodiode array rotated 
360° around the incident light was used to detect the total 
intensity of the reflected and scattered light. The results 
indicated that the peak and the valley of the periodic inten¬ 
sity distribution corresponded to the rotation angle perpen- 
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dicular and parallel to the grain direction, respectively. “ 
Although it is an effective technology, the method proposed 
is expensive and time-consuming. To improve this method, 
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Fig. 1. Schematic diagram of the 
experimental system 
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the laser light scattering pattern was used to determine the 
grain angle and diving angle. The light scattering pattern 
image was a quasi-ellipse caused by the tracheid effect. 
Principal component analysis (PCA) was used to determine 
the grain angle and the ratio of the minor-to-major axes of 
the ellipse was used to evaluate the diving angle. 8 

In our research, a new ellipse analysis of the laser light 
scattering pattern image was proposed to accurately deter¬ 
mine the ellipse parameters automatically and thereby de¬ 
termine the grain angle. 


Materials and methods 

The schematic diagram of the system is shown in Fig. 1. 
It consists of a He-Ne laser source (Edmund Optics, 
NT61332), a 10X Beam Expander (Edmund Optics, NT 
55578), a cubic beam splitter (Edmund Optics, NT 32503), a 
CCD-camera (Victor KY-F350, Victor JVC), a rotatable 
table, some ancillary holding parts, and a host computer 
fitted with an image processing program developed and 
installed by the authors. The materials selected for the ex¬ 
periment were sugi (Cryptomeria japonica D.Don) and 
Japanese beech (Fagus crenata blume). Samples were 
placed on the rotatable table. The laser light propagating 
path from the laser source to the CCD-camera is shown in 
Fig. 1. Considering the working distance and processing 
efficiency, the scattering pattern image was captured at 
150 dpi (dot per inch) by adjusting the working distance and 
the focus length of the CCD camera. The computer pro¬ 
gram that was used to process the laser light scattering 
pattern image and to analyze it by means of modified 
Hough transform ellipse analysis was written in Visual 
Basic 6.0 as code resource. The flow chart is shown in Fig. 2. 



Image processing algorithm 

After capturing and inputting the laser light scattering im¬ 
age, the Sobel edge detecting operation 9 was performed to 


Fig. 2. Flow chart of the computer program tasks 
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detect the elliptical edge points. This operation can be de¬ 
scribed by the following formulas: 
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where k is the iterating count. After definite iteration times, 
Xf +1) will converge quickly to its accurate solution with any 
arbitrary initials Xf ] (j = 1, 2,..., 5). 

For n elliptic edge points, all of the combinations of the 
five-point pair can be denoted as (p b p p p k , p b p m ) (i ^j^k 
^ ^ m, 1 < i, j, k, /, m < n). However, fitting all combina¬ 
tions of five points will be overly time-consuming. To save 
computing time, the MHT operation is performed on all 
regularly spaced five-point pairs denoted as follows: 


where S v and S v are the convolution masks of x and y direc¬ 
tional intensity gradients; A is the 3X3 neighborhood 
matrix with the geometric center of the transformed point 
(x, y); S xy is the desired pixel value of the point (x, y). If the 
value of S xy is greater than a threshold determined by pre¬ 
liminary experiments, the point (x, y) will be detected as 
one of the elliptical edge points. 

In order to find the four ellipse parameters (major axis, 
minor axis, orientation, and eccentricity), which are equiva¬ 
lent to determining the grain angle, all the elliptical edge 
points were passed through the ellipse analysis of the 
modified Hough Transform (MHT). 1012 The adopted MHT 
algorithm can be described as follows. 

A general representation of an ellipse is denoted by 
Eq. 4. 

Ax 2 + Bxy + Cy 2 + Dx + Ey + 1 = 0 (4) 

with the constraint of 4 AC — B > 0. An ellipse is uniquely 
defined by a minimum of five points. The conic coefficients 
of A, B, C, D , E can be determined by substituting the five 
points into Eq. 4 and solving the five simultaneous equa¬ 
tions shown as follows. 

Ax? + Bx l y l + Cy 2 + Dx t + Ey t + 1 = 0 (5) 

where x h y t (i = 1,2,..., 5) is any arbitrary point’s coordi¬ 
nates of the elliptic edge detected by the above Sobel opera¬ 
tion. Note that all of x h y / (/ = 1, 2,..., 5) are becoming 
coefficients while all of the ellipse coefficients (A, B, C, D, 
E) are becoming variables in the above linear simultaneous 
equations. To be easily understood, the linear simultaneous 
equations can be rewritten as Eq. 6. 

5 

Y.a.Xj = -1, i = 1,2,... ,5 (6) 

7=1 


(T/5 Pi+s') Pi+2si Pi+3si Pi+As ) 

= 1,2 ,...,n — 4s; s = 1,2,..., (n — l)/4^. 



s is the separation between points. Note that the maximum 
computational complexity is being reduced. Each ellipse 
fitting through a five-point pair (p b p i+s , p i+2s , p i+3s , p i+4s ) 
provides an estimate (A is , B is , C is , D is , E is ) of the conic 
coefficients. For some sets of five points, the conic may be a 
parabola or hyperbola instead of an ellipse. If the deduced 
conic coefficients cannot satisfy the additional constraint of 
4-4; > 0, these five nonelliptic points will be re¬ 

jected. If the conic coefficients can satisfy the constraint of 
4/l,.C'„ - /»•'. > 0, these coefficients, denoted as (A( tS , B ’ Us , 


C'i,s, A Is, E ls) to differ from (A, s , B is , C, s , D is , EJ, will be 
used to compute the desired four ellipse parameters de¬ 
scribed by the following formulas: 
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( 11 ) 


e itS = aus/Ks (I 2 ) 

where 0 is is the ellipse orientation; a is , b is are the semi 
minor and major axes of the ellipse; and e is is the ellipse 
eccentricity. The intensity histograms of the above ellipse 
parameters were computed. The final four parameters were 
taken as the dominant peaks of each set of parameters 


where a tj are coefficients of the linear equations; Xj are 
variables of the linear simultaneous equations; i is the total 
number of variables in the linear simultaneous equations; j 
is the total number of linear simultaneous equations. 

The linear simultaneous equations are solved by the 
Gauss-Seidel iteration algorithm 13 described as follows: 
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where 6 is the final accumulating orientation of the ellipse, 

y\ 

which can be used to determine the grain angle; a and b are 
the final accumulating semimajor and/or minor axes of the 
ellipse, respectively; and e is the final accumulating eccen¬ 
tricity of the ellipse that can be used to determine the diving 
angle. 8 


Results and discussion 

The picture of the laser light scattering spot is shown in 
Fig. 3. These pictures were captured at every rotating posi¬ 
tion for the same location on the sample surface while the 
planed sugi sample was rotated every 15° from 0° to 360°. 
The horizontal direction was defined as the base axis for 
determining the grain angle and the horizontal center line 
of the 0° picture was defined as the base rotating axis. As 
shown in Fig. 3, all of the laser spot scattering patterns are 
quasi-ellipsoidal and their orientations changed with the 
rotation. This phenomenon is considered to be caused by 
the tracheid effect. To prove this, base pictures were also 
captured. A picture of the sample surface was taken using a 
color CCD camera and the picture was printed onto a high 


quality white paper by a color printer. The picture of the 
laser spot scattering pattern, which was called the base pic¬ 
ture, was captured at the same position on the printed paper 
as the real sample was operated. It can be readily seen that 
the base picture is a quasi-circle; the effect of the tracheid is 
manifested well to give a quasi-ellipse in the captured pic¬ 
ture of the laser spot scattering pattern. 

The quasi-ellipse edge of the picture shown in Fig. 3 
detected by the Sobel operation is shown in Fig. 4. All of the 
edge points were detected correctly. Although the laser 
spot scattering pattern is not a clear-out ellipse, most of its 
edge points tend to share the same ellipse parameters. The 
changing orientation of the laser spot scattering pattern 
with the rotation is clearer than that in Fig. 3. The quasi¬ 
ellipse edge points were analyzed by the ellipse analysis 
technology of the modified Hough transform. The param¬ 
eters of major and minor axes are shown in Fig. 5 and the 
eccentricity parameter is shown in Fig. 6. As shown in Figs. 
5 and 6, all the parameters are approximately the same for 
each measuring angle. According to the planed sugi sample, 
the major axis fluctuates around 86 pixels (or about 15 mm) 
with maximum deviation of ±1.4 pixels (or about 0.2 mm) 
and the minor axis fluctuates around 49 pixels (or about 
8 mm) with maximum deviation of ±1.2 pixels (or about 


Fig. 3. The laser spot scattering 
pattern pictures. Rotation angles 
are stated relative to the horizon¬ 
tal plane 
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Fig. 4. The detected quasi-ellipses 
of the laser spot scattering pattern. 
Each binary image corresponds to 
the image in the same positions in 
Fig. 3 




Measuring position (°) Fig. 6. Eccentricities of the ellipses 

Fig. 5. The major and minor axes of the ellipses. Symbols and vertical 
bars denote the mean values and standard deviations, respectively 
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Fig. 7. The grain angle measured automatically 


0.2mm); the eccentricity values fluctuate around 0.58. The 
results, which show that all the quasi-ellipses share similar 
ellipse parameters, indicate that the adopted ellipse analysis 
technology of the modified Hough transform is efficient. 
The reason is that all of the edge points are treated as points 
of ellipse sets instead of an exact ellipse and all five ellipse 
parameters are computed from the parameter population 
distributions of those ellipse sets. 

The grain angle measured automatically by the method 
proposed in this research is shown in Fig. 7. The measured 
grain angle coincided with the real grain angle and there is 
an almost linear relation between the measured grain angle 
and the real grain angle. Some researchers have reported 
that the wood species and surface roughness affected the 
measuring accuracy of grain angles using optical scatter¬ 
ing. 14 To study the effect of the surface quality on the effi¬ 
ciency of the proposed method in our research, Japanese 
beech and sugi were used with rough-sawn and planed fin¬ 
ishes. There was no significant difference between the mea¬ 
sured grain angles of Japanese beech and those of sugi 
under the above two machining conditions. The deviation 
of the measured grain angle in the planed samples was less 
than that in rough-sawn samples for both Japanese beech 
and sugi. The deviations for sawn samples were 4.6° and 
2.9° for sugi and Japanese beech, respectively; the devia¬ 
tions for planed samples were 0.46° and 1.3° for sugi and 
Japanese beech, respectively. The results show that the sur¬ 
face quality affects the accuracy of the system. However, 
the effect on the final accuracy is not significant. In addition, 
the measured grain angle of sugi was more accurate than 
that of Japanese beech for planed samples and the mea¬ 
sured grain angle of Japanese beech was more accurate than 
that of sugi for sawn samples. This is because the interaction 
of the different tracheids and the different machining condi¬ 
tions causes a different blurred contour of the laser spot. 
The difference between them is not significant and the ef¬ 
fects of the interaction between the different sample prop¬ 
erties and the different machining conditions are generated 
only during the process of generating pictures. After gener¬ 
ating and capturing the laser spot scattering pattern picture 


accurately, the accuracy of the whole system depends on the 
image processing technology. 


Conclusions 

We developed an on-line system to determine the grain 
angle using the tracheid effect to generate the laser scatter¬ 
ing pattern image, which was analyzed by the new ellipse 
analysis of the modified Hough transform. The results have 
shown that the on-line image analyzing system is effective 
and sufficiently robust to determine the grain angle for 
both sugi and Japanese beech for rough-sawn and planed 
samples. The measured grain angle is proportional to the 
real grain angle with zero intercept and a slope of 1. The 
deviations of sugi were 4.6° and 0.46° for sawn and planed 
samples, respectively; the deviations of Japanese beech 
were 2.9° and 1.3° for sawn and planed samples, respec¬ 
tively. There is an insignificant difference between the re¬ 
sults for Japanese beech and sugi under two machining 
conditions. However, the measured grain angle of sugi was 
better than that of Japanese beech for planed samples and 
the measured grain angle of Japanese beech was better than 
that of sugi for sawn samples. The results of both sugi and 
Japanese beech for planed samples were better than those 
for sawn samples. 
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